Nucleation of A^Zr and AI3SC in aluminum alloys: 
from kinetic Monte Carlo simulations to classical theory. 
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Zr and Sc precipitate in aluminum alloys to form the compounds AbjZr and AI3SC which for low 
supersaturations of the solid solution have the LI2 structure. The aim of the present study is to 
model at an atomic scale this kinetics of precipitation and to build a mesoscopic model based on 
classical nucleation theory so as to extend the field of supersaturations and annealing times that 
can be simulated. We use some ab-initio calculations and experimental data to fit an Ising model 
describing thermodynamics of the Al-Zr and Al-Sc systems. Kinetic behavior is described by means 
of an atom- vacancy exchange mechanism. This allows us to simulate with a kinetic Monte Carlo 
algorithm kinetics of precipitation of AbjZr and AI3SC. These kinetics are then used to test the 
classical nucleation theory. In this purpose, we deduce from our atomic model an isotropic interface 
free energy which is consistent with the one deduced from experimental kinetics and a nucleation 
free energy. We test different mean-field approximations (Bragg- Williams approximation as well as 
Cluster Variation Method) for these parameters. The classical nucleation theory is coherent with the 
kinetic Monte Carlo simulations only when CVM is used: it manages to reproduce the cluster size 
distribution in the metastable solid solution and its evolution as well as the steady-state nucleation 
rate. We also find that the capillary approximation used in the classical nucleation theory works 
surprisingly well when compared to a direct calculation of the free energy of formation for small LI2 
clusters. 
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I. INTRODUCTION 

Precipitation kinetics of a metastable solid solution is 
known to be divided in three successive stages: the nu- 
cleation, growth, and coarsening of nuclei of the new sta- 
ble phase. The first stage of precipitation is of great 
practical interest but difficult to observe experimentally. 
Kinetic Monte Carlo simulation is the suitable tool for 
a numerical prediction of a nucleation kinetics**^ but a 
rationalization of the results is difficult and atomic sim- 
ulations cannot reach very low supersaturations. On 
the other hand, classical descriptions of these different 
stages*^ are well established and the associated mod- 
els are now widely used to understand experimental ki- 
netics and to model technological processe a 5 i 6 i 7 ' 8 . Re- 
cently, classical nucleation theory has been shown to be 
in good agreement with more reliable atomic models by 
way of a direct comparison with kinetic Monte Carlo 
simulations***i**i***. These different studies included 
decomposition of a metastable solid solution for a demix- 
ing binary system on a surface* or in the bulki*** and 
kinetics of electrodeposition on a surfacei*. In this last 
study, Berthier et al. show that physical parameters of 
classical nucleation theory have to be carefully calculated 
so as to reproduce atomic simulations. In the present 
paper, we want to extend the range of comparison be- 
tween classical nucleation theory and atomic simulations 
by studying the case of an ordering system on a frustrated 
lattice. We thus choose to model kinetics of precipitation 
of a LI2 ordered compound formed from a solid solution 
lying on a face centered cubic (fee) lattice. 

For fee lattice it is now well established that one has 
to use a mean field approximation more accurate than 



the widely used Bragg- Williams one in order to calcu- 
late thermodynamic properties**. The Cluster Variation 
Method (CVM)i*>i* enables one to obtain phase dia- 
grams which are in quantitative agreement with thermo- 
dynamic Monte Carlo simulations**!*. When CVM is 
used, frustration effects on the tetrahedron of first near- 
est neighbors and short range order due to interactions 
are considered in a satisfying way enabling one to pre- 
dict quantitatively thermodynamic behavior. Neverthe- 
less, the use of CVM is often restricted to the calculation 
of equilibrium properties and, thus, for computing ther- 
modynamic properties of the metastable supersaturated 
solid solution in classical nucleation theory one merely 
considers Bragg- Williams approximation. The purpose 
of this paper is then to show that the use of CVM calcula- 
tions with classical nucleation theory leads to a satisfying 
description of the metastable solid solution and extend 
the range of supersaturations that can be modelled with 
this theory. 

In this purpose we build an atomic model which al- 
lows us to study kinetics of precipitation of A^Zr and 
AI3SC. The two considered binary systems, Al-Zr and 
Al-Sc, have different kinetic properties: the interaction 
with vacancies is repulsive for Zr atoms whereas it is at- 
tractive for Sc atoms. On the other hand, for low super- 
saturations, thermodynamics of both systems are quite 
similar. A^Zr has the stable DO23 structure!*, but for 
small supersaturations of the solid solution, AI3Z1- precip- 
itates with the metastable LI2 structure and precipitates 
with the DO23 structure only appear for prolonged heat 
treatment and high enough supersaturations****!. As 
for AI3SC, the stable structure is LI^A* and thus only 
LI2 precipitates have been observed during experimen- 
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tal kinetics^SiSiSi. In this study we mainly focus on the 
nucleation stage and therefore we consider that both Zr 
and Sc lead to the precipitation of a compound having 
the LI2 structure. In this context, Al-Zr and Al-Sc sys- 
tems are really similar from a thermodynamic point of 
view unlike their kinetic behavior. It is then interest- 
ing to study these two systems in parallel and to see if 
classical nucleation theory manages to reproduce atomic 
simulations for these two different kinetic behaviors. 

The atomic model used in kinetic Monte Carlo simula- 
tions is built using experimental data as well as ab-initio 
calculations. We deduce from it physical parameters en- 
tering mesoscopic models like classical nucleation theory 
and show how this theory compares to atomic simula- 
tions for different supersaturations and different anneal- 
ing temperatures. The capillary approximation used in 
classical nucleation theory is then discussed as well as dif- 
ferent mean field approximations that can be combined 
with it. 



II. ATOMIC MODEL 



TABLE I: First and second nearest neighbor pair effective en- 
ergies (in eV). Only interactions different from zero are pre- 
sented. 



A1A1 

e (1) 

fc ZrZr 

e (1) 

c ScSc 
fc AlZr 
fc AlSc 

e (1) 
E vv 

(2) 
fc AlZr 

J2) 



-0.560 
-1.045 
-0.650 

-0.979 + 24.4 x 10" 6 T 
-0.759 + 21.0 x 10 _6 T 
-0.084 

+0.101 - 22.3 x 10" 6 T 
+0.113 - 33.4 x 10" 6 T 



e (1) 

e AlV 
fc ZrV 



-0.222 
-0.350 
-0.757 



With such a model, as long as vacancy concentration 
can be neglected, thermodynamic behavior of Al-X sys- 
tem (X = Zr or Sc) only depends on the order energies 
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A. Al-Zr and Al-Sc thermodynamics 

In order to simulate thermodynamic behavior of Al-Zr 
and Al-Sc binary systems, we use a rigid lattice: config- 
urations of the system are described by the occupation 
numbers p l n with p\ = 1 if the site n is occupied by an 
atom of type i and p % n = if not. Energies of such con- 
figurations are given by an Ising model with first and 
second nearest neighbor interactions. This is the sim- 
plest model to simulate precipitation of a stoichiometric 
AI3X compound in the LI2 structure. Indeed one has to 
include second nearest neighbor interactions, otherwise 
LI2 precipitates do not show perfect AI3X composition. 
On the other hand, there is no use to consider interac- 
tions beyond second nearest neighbors as these interac- 
tions are significantly lower than first and second nearest 
neighbor interactions^ 5 -. We could have considered inter- 
actions for clusters other than pairs too, but we showed 
that the use of interactions for first nearest neighbor tri- 
angle and tetrahedron does not really change the kinetics 
of precipitation^: the Onsager coefficients defining diffu- 
sion in the solid solution are unchanged with or without 
these interactions as well as the nucleation free energy. 
Thus, in our model, the energy per site of a given config- 
uration is 



We use experimental data combined with ab-initio calcu- 
lations to obtain these order energies for Al-Zr and Al-Sc 
systems. 

First nearest neighbor order energies are chosen so as 
to correctly reproduce formation energies of AbjZr and 
Al 3 Sc compounds in Ll 2 structure, AF(A1 3 X, LI 2 ) = 
3wW. For Al 3 Zr, we use the free energy of formation 
that we previously calculated 2 ^. For AI3SC, we calculate 
the enthalpy of formation with the full-potential linear- 
muffin-tin-orbital method 2 in the generalized grad ient 
approximation 2 ^ and we use the value of Ref. |29| and 
l30l for the vibrational contribution to the free energy of 
formation: 



AF(Al 3 Zr, Ll 2 ) 
AF(Al 3 Sc, Ll 2 ) 



= -0.530 + 73.2 x 10~ 6 T eV, 
= -0.463 + 62.9 x I0~ 6 T eV. 



Second nearest neighbor interactions are chosen so as 
to reproduce Zr and Sc solubility limits in Al. Indeed 
these limits only depend on order energy u/ 2 \ as can be 
seen from the low temperature expansion^! to the second 
order in excitation energies: 



exp 



(-6w (2) /fcT) +6exp(-10w (2) /fcT). (4) 
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where the first and second sums respectively runs on all 
first and second nearest neighbor pairs of sites, N s is the 



number of lattice sites, and e\"-' are the respective 
effective energies of a first and second nearest neighbor 
pair in the configuration 



,( 2 ) 



We check using the CVM in the tetrahedron-octahedron 



approximatio; 



.15.16 



that this low temperature expansion 



for the solubility limit is correct in the whole range of 
temperature of interest, i.e. until Al melting tempera- 
ture (T mel — 934K). For Al-Zr interactions, as we want 
to model precipitation of the metastable LI2 structure of 
Al 3 Zr compound, we use the metastable solubility limit 
that we previously obtained from ab-initio calculations 25 , 
whereas for Al 3 Sc the LI2 structure is stable and we use 



3 



the solubility limit arising from a thermodynamic mod- 
elling of experimental data^: 

x% = exp ((-0.620+ 155 x 10~ 6 T) eV/kT), 
x e s l = exp ((-0.701 + 230 x 10~ 6 T) eV/kT). 

One should notice that these solubility limits have been 
found to be consistent with ab-initio calculations^^ and 
thus with the formation energies we used for A^Zr and 
and Al 3 Sc. 

Unlike thermodynamics, kinetics does not only de- 
pend on order energies but also on effective energies 
and . We deduce them from w' 1 ' and by 
using experimental values for cohesive energies of pure 
elements 3 ^: E coh (A\) = 3.36 eV, E coh (Zr) = 6.27 eV, 
and I? (Sc) — 3.90 cV. We assume that second near- 
est neighbor interactions do not contribute to these co- 

(2) 

hesive energies, i.e. e^x = (X = Al, Zr, or Sc) and 
we neglect any possible temperature dependence of these 
energies. Therefore, the cohesive energy of X element is 
E coh (X) = 6exx- Resulting effective energies are pre- 
sented in table CI 



B. Al-Zr and Al-Sc kinetics 

We introduce in the Ising model atom- vacancy interac- 
tions for first nearest neighbors (Tab.QJ, so as to consider 
the electronic relaxations around the vacancy. Without 
these interactions, the vacancy formation energy E v ° r in 
a pure metal would necessarily equal the cohesive energy 
which is in contradiction with experimental data, e^iv 
and £zrv are deduced from vacancy formation energy re- 
spectively in pure Al24, Ey OT — 0.69 eV, and in pure Zv^, 
E y or = 2.07 eV. For Zr, this energy corresponds to the 
hep structure which is quite similar to the fee one (same 
first nearest neighborhood). Therefore, we assume that 
the vacancy energy is the same in both structures. This 
is then possible to correct this formation energy to take 
into account the difference between Al and Zr equilibrium 
volumes, but this leads to a correction of ~ 10% for E v ° r 
and does not really change the physical interaction be- 
tween Zr atoms and vacancies. We thus choose to neglect 
such a correction. To compute the interaction between 
Sc atoms and vacancies, we can directly deduce CgJv 
from the experimental binding energy in aluminum 3 ^, 

#scv = 4iv + - 4cv - 4iai = 0-35 eV at 650 K. 
Such an experimental data does not exist for Zr impurity, 
but we can check that the physical interaction we obtain 
is correct. The binding energy deduced from our set of 
parameters is strongly negative (S|"y = —0.276 eV at 
650 K). This is in agreement with the experimental fact 
that no attraction has been observed between vacancy 
and Zr impurity 3 ^!. This repulsion in the case of Zr im- 
purity and this attraction in the case of Sc impurity are 
related to the difference of cohesive energies between Zr 



TABLE II: Kinetic parameters: contribution of the jumping 
atom to the saddle point energy e a J and attempt frequency 
u a for a = Al, Zr, and Sc atoms. 



e Al — 


-8.219 eV 


"Al = 


1.36 x 10 14 


Hz 


P S P - 
e Zr — 


-11.464 eV 


^Zr - 


9 x 10 16 


Hz 


C Sc — 


-9.434 eV 


^Sc = 


4 x 10 15 


Hz 



and Sc, showing thus that elastic relaxations around the 
vacancy are not the dominant effect. It could explain why 
Zr diffusion coefficient in aluminum is so low compared 
to the Sc one. Some ab-initio calculations have been 
made to compute this binding energy with a vacancy for 
all transition metals in aluminum 38 . They obtained in 
the case of Zr as well as Sc impurity a repulsive inter- 
action with a vacancy. This is in contradiction with the 
experimental data we use to compute Sc-vacancy inter- 
action. Such a disagreement may arise from approxima- 
tions made in the calculation (Kohn-Korringa-Rostoker 
Green's function method) as the neglect of atom relax- 
ations and the box that includes only the first nearest 
neighbors of the impurity- vacancy complex. Nevertheless 
these ab-initio calculations showed that a binding energy 
as large as 0.35 eV is possible as the value obtained for 
Sr impurity was even larger. 

We use the experimental value of the divacancy binding 
energy 3 ^, E^v — 0.2 eV, in order to compute a vacancy- 
vacancy interaction, ey V = 2e A 1 1 y — e^iAl — ^2V ■ ^ we 
do not include this interaction and set it equal to zero 
instead, we obtain a binding energy which is slightly too 
low, divacancies being thus not as stable as they should 
be. Some recent ab-initio calculations 3 ^^ have shown 
that divacancies should be actually unstable, the non- 
Arrhcnius temperature dependence of the vacancy con- 
centration arising from anharmonic atomic vibrations. 
Nevertheless, this does not affect our Monte Carlo simu- 
lations as we only include one vacancy in the simulation 
box, but this divacancy binding energy should be con- 
sidered more seriously if one wants to build a mean field 
approximation of our diffusion model or if one wants to 
compensate vacancy trapping by adding new vacancies 
in the simulation box. 

Diffusion is described through vacancy jumps. The 
vacancy exchange frequency with one of its twelve first 
nearest neighbors of type a is given by 

T Q -v = v a exp ( — J , (5) 

where v a is an attempt frequency and the activation en- 
ergy E® ct is the energy change required to move the a 
atom from its initial stable position to the saddle point 
position. It is computed as the difference between the 
contribution e^f of the jumping atom to the saddle point 
energy and the contributions of the vacancy and of the 
jumping atom to the initial energy of the stable posi- 
tion. This last contribution is obtained by considering 
all bonds which are broken by the jump. 
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The attempt frequency v a and the contribution e* p of 
the jumping atom to the saddle point energy can depend 
on the configuration 4 ^^ 2 * 4 ^. Nevertheless, we do not have 
enough information to see if such a dependence holds 
in the case of Al-Zr or Al-Sc alloys. We thus assume 
that these parameters depend only on the nature of the 
jumping atom. We fit the six resulting kinetic parameters 
(Tab.[nJ) so as to reproduce Al self-diffusion coefficient 44 
and Zr^ 4 ^ and So^& impurity diffusion coefficients: 

Dai* = 0.173 x l(T 4 exp (-1.30 eV/kT) m 2 .s _x , 
£>zr* = 728 x l(T 4 exp (-2.51 eV/fcT) m 2 .s _1 , 
D Sc . = 5.31 x lCT 4 exp (-1.79 eV/fcT) m 2 .s _1 . 

Zr is diffusing slower than Sc which itself is diffusing 
slower than Al. The difference between diffusion coef- 
ficients is decreasing with the temperature, but at the 
maximal temperature we considered, i.e. T = 873 K, it 
still remains important as we then have -Dai* ~ 20Z?s c * ~ 
200dD Zr *- 



C. Monte Carlo simulations 

We use residence time algorithm to run kinetic Monte 
Carlo simulations. The simulation boxes contain N s = 
100 3 or 200 3 lattice sites and a vacancy occupies one of 
these sites. At each step, the vacancy can jump with one 
of its twelve first nearest neighbors, the probability of 
each jump being given by Eq. lf5|. The time increment 
corresponding to this event is 

1 1 
'-^(1-13^)CV(A1) E^rlVv' (6) 

where is the nominal concentration of the simulation 
box (X = Zr or Sc) and CV(A1) the real vacancy con- 
centration in pure Al as deduced from energy parame- 
ters of table [I] The first factor appearing in Eq. © is 
due to the difference between the experimental vacancy 
concentration in pure Al and the one observed during 
the simulations. The dependence of this factor with the 
concentration i;| reflects that for each impurity the cor- 
responding lattice site and its twelve first nearest neigh- 
bors cannot be considered as being pure Al. It is correct 
only for a random solid solution in the dilute limit, but 
concentrations considered in this study are low enough 
so the same expression can be kept for this factor. The 
absolute time scale is then obtained by summing only 
configurations where the vacancy is surrounded by Al 
atoms in its first nearest neighborhood, i.e. where the 
vacancy is in pure Al. This ensures that the influence on 
the time scale of the thermodynamic interaction between 
the vacancy and the Sc or Zr impurity is correctly taken 
into account. The method employed here is equivalent to 
measuring the fraction of time spent by the vacancy in 
pure Al during the simulation and multiplying then all 
time increments by this factor as done in Ref. ^3- We 



check by running Monte Carlo simulations for different 
sizes of the box that results do not depend on the effective 
vacancy concentration. 

For low impurity concentration, residence time algo- 
rithm can be sped up by noticing that in most of the ex- 
plored configurations the vacancy is located in pure Al, 
i.e. on a lattice site where all exchange frequencies are 
equal to the one in pure Al. In such a configuration, the 
move of the vacancy can be associated to a random walk 
and the corresponding time increment is known a priori. 
Lattice sites corresponding to such configurations can be 
detected in the beginning of the simulation and the cor- 
responding tables need to be modified only each time 
the vacancy exchanges with an impurity 4 ^. For impurity 
concentration in the range 5 x 10~ 3 < < 1 x 10~ 2 , 
the algorithm is sped up by a factor ~ 2. This allows us 
to simulate lower supersaturations of the solid solution 
than we could have done with a conventional algorithm. 

So as to follow kinetics of precipitation in the simula- 
tion box, we have to give us a criterion to discriminate 
atoms belonging to the solid solution from those in LI2 
precipitates. As, according to the phase diagram, the 
stoichiometry of theses precipitates is almost perfect, we 
only look at Zr or Sc atoms and assume for each of these 
atoms in a LI 2 cluster that three associated Al atoms 
belong to the same cluster. Zr (or Sc) atoms are counted 
as belonging to a cluster having LI2 structure if all their 
twelve first nearest neighbors are Al atoms and at least 
one of their six second nearest neighbors is a Zr (or Sc) 
atom. This criterion works only for dimers and bigger 
clusters and then all Zr or Sc atoms not belonging to such 
clusters are considered to be monomers. We only counted 
as precipitates LI2 clusters bigger than a critical size, 
i.e. containing more Zr or Sc atoms than a critical num- 
ber n^, this critical number being chosen as the initial 
one given by classical nucleation theory (cf. section ITTTll . 
Clusters smaller than this critical size are unstable and 
will re-dissolve into the solid solution. Therefore atoms 
contained in such clusters are counted as belonging to 
the solid solution. With this criterion, we can measure 
during the atomic simulations the number of stable pre- 
cipitates and their mean size (Fig. and [2J. The solid 
solution concentration is then defined at each step by the 
relation 

xx = ^ n xC"„ x , (?) 

«x = l 

where C„ x is the number of LI2 clusters containing nx Zr 
or Sc atoms normalized by the number of lattice sites, i.e. 
the instantaneous probability to observe such a cluster in 
the simulation box. 

All starting configurations for simulations are com- 
pletely disordered (random) solid solutions. We thus sim- 
ulate infinitely fast quenching from high temperatures. 
During the first steps of the precipitation, the number 
of stable precipitates is varying quite linearly with time 
(Fig. |2l • The slope of this linear relation gives a measure 



5 




FIG. 1: Monte Carlo simulation of the kinetics of precipitation of AI3SC for a supersaturated aluminum solid solution of nominal 
concentration x$ c — 0.005 at T = 773 K. The simulation box contains 8 x 10 6 lattice sites. Only Sc atoms belonging to LI2 
precipitates are shown. The critical size used is ng c = 13. 




f 

-t 1 1 1 1 1 

0.1 0.2 0.3 

f(s) 



FIG. 2: Kinetics of precipitation of a supersaturated alu- 
minum solid solution of nominal concentration Kg c = 0.005 
at T = 773 K: evolution with time of the number N 3p of 
stable precipitates in the simulation box (normalized by the 
number of lattice sites N s ), of stable precipitate average size 
< nsc >s P , and of Sc concentration xsc in the solid solution. 
The critical size used to discriminate stable precipitates from 
sub-critical clusters is n Sc = 13. Some of the corresponding 
simulation configurations are shown in Fig. 



of the steady-state nucleation rate J st , i.e. the number 
of stable precipitates appearing by time unit during the 
nucleation stage. 



III. CLASSICAL NUCLEATION THEORY 

In order to compare kinetic Monte Carlo simulations 
with classical nucleation theory, we need to define the 
formation free energy AG n (x x ) of a LI2 cluster contain- 
ing n atoms (n = 4nz r or 4ng c ) embedded in a solid 
solution of nominal concentration x^. Usually, one uses 
the capillary approximation and considers a volume con- 
tribution, the nucleation free energy AG nuc (x x ), and a 
surface contribution corresponding to the energy cost to 
create an interface between the solid solution and the LI2 
precipitate, 

/q X 1 / 3 

AG„(4) = nAG nuc (x° x ) + ( ^ J n 2 / 3 a 2 a, (8) 

where a is the lattice parameter and a the interface free 
energy. For a supersaturated solution, AG„(ix) shows a 
maximum in n* = 4n| r or 4rig c corresponding to the crit- 
ical size used to follow the kinetics of precipitation during 
the Monte Carlo simulations (c/. section III Cll . We now 
have to calculate the nucleation free energy AG nuc (x x ) 
and the interface free energy a corresponding to the set 
of atomic parameters presented in the previous section. 

A. Nucleation free energy 

The nucleation free energy to precipitate AI3X (X = 
Zr or Sc) is^ 

AG— (4) = ^(mai(4 9 )-Mai(4)) 

+i( Mx (40-^x(4)), (9) 

where /xai(^x) an d fJ-x(%x) are the chemical potentials of 
respectively Al and X components in the solid solution 
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FIG. 3: Dependence with temperature of the free energies of 
the solid solution / AbjZr (top) and solid solution / AI3SC 
(bottom) interfaces for the [100], [110], and [111] directions, 
and associated isotropic free energy a obtained from Wulff 
construction. 



of concentration xx, is the equilibrium concentration 
of the solid solution, and x x the nominal concentration. 
The factors 3/4 and 1/4 arise from the stoichiometry of 
the precipitating phase AI3X. We use the CVM in the 
tetrahedron-octahedron approximationi^^ to calculate 
chemical potentials entering Eq. @. This is the min- 
imum CVM approximation that can be used with first 
and second nearest neighbor interactions. Within this 
approximation, all correlations inside the tetrahedron of 
first nearest neighbors and the octahedron linking the 
centers of the six cubic faces are included in the calcu- 
lation of the chemical potentials. Usually one does not 
consider these correlations in the calculation of the nu- 
cleation free energy and merely uses the Bragg- Williams 
approximation to obtain AG nuc but we will see in section 
IIV Bl that this leads to discrepancies between results of 
atomic simulations and predictions of classical nucleation 
theory. 



B. Interface free energy 

1. Plane interfaces 

We calculate interface free energies between the alu- 
minum solid solution assumed to be at equilibrium close 
to the interface and the LI2 precipitates for three dif- 
ferent directions of the interface ([100], [110], and [HI]). 
If phases are assumed to be pure, the different interface 



energies are simply related by the equation 

1 1 ujM , . 

a wo = —oixo = ^111 = — • (10) 

At finite temperature, one has to consider that the solid 
solution is not pure Al and that the LI2 structure differs 
from AI3X stoichiometry. Moreover, to minimize the en- 
ergy cost due to the interface, concentrations and order 
parameters of planes near the interface can differ from 
those in the bulk. So as to take into account such a re- 
laxation, we calculate these interface free energies within 
the Bragg- Williams approximation. A better statistical 
approximation based on CVM is too cumbersome and we 
only check for the [100] direction that we obtain the same 
value of the free energy in the whole range of tempera- 
tures with a CVM calculation in the tetrahedron approx- 
imation. 

At finite temperature, we still observe that cioo < 
ciio < fni (Fig|3J). Nevertheless, as the relaxation is 
small for the interface in the [100] direction and impor- 
tant for the [111] direction, the difference between in- 
terface energies is decreasing with temperature. This 
indicates that precipitates are becoming more isotropic 
at higher temperatures. Using Wulff constructionist to 
determine the precipitate equilibrium shape, we find that 
precipitates will mainly show facets in the [100] direc- 
tions and that facets in the [110] and [111] directions are 
small but becoming more important with increasing tem- 
perature. Comparing these predicted equilibrium shape 
with the ones observed during the atomic simulations, 
we find a good agreement (Fig. at low temperatures 
(T ~ 723 K), precipitates are cubic with sharp [100] in- 
terfaces, whereas at higher temperatures interfaces are 
not so sharp. For AI3SC, Marquis and SeidmanS ex- 
perimentally observed precipitates showing facets in the 
[100], [110], and [111] directions at T = 573 K, with [100] 
facets tending to disappear at high temperatures. This 
is well reproduced by our atomic model, the main differ- 
ence being that the experimentally observed [100] facets 
are less important compared to the other ones than in 
our study. 

Asta et til Ml used a cluster expansion of ab-initio calcu- 
lations to obtain the same interface energies in Al-Sc sys- 
tem. The energies they got are higher than ours: a 2 trioo 
is varying from 167 to 157 meV between OK and the melt- 
ing temperature (T mel = 934K) and a 2 cr m from 233 to 
178 meV. The difference could be due to the limited range 
of our interactions compared to Asta's ones. This could 
explain too why the interface free energies we obtain are 
decreasing more rapidly with temperature especially in 
the [111] directional. 

Hyland et alJH also calculated interface free energies 
with an empirical potential for Al-Sc, but the energies 
they obtained are really low compared to ours and Asta's 
ones as well as compared to the isotropic interface free 
energy they measured^. Some of the discrepancy can be 
due to relaxations of atomic positions which are consid- 
ered in their study and are missing in our. It may also 
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(a) T = 723 K 




(b) T = 873 K 



FIG. 4: AI3SC precipitate observed during Monte Carlo sim- 
ulations at different temperatures and corresponding Wulff 
construction obtained from the interface free energies calcu- 
lated at the same temperatures. For Monte Carlo simula- 
tions, only Sc atoms are shown. Atom color corresponds to 
the number of Sc atoms as second nearest neighbors. For a 
[100] interface, it should be 5, for a [110] 4, and for a [111] 3. 



indicate that the potential they used is not really well 
suited to describe solid solution / AhjSc interfaces. 



2. Average interface free energy 

We use the Wulff constructionist to define an 
isotropic free energy a from the free energies ctioo, ono, 
and <7ioo- v is defined so as to give the same interface 
free energy for a spherical precipitate having the same 
volume as the real faceted one. Details of calculations 
can be found in appendix B The free energy a is higher 

than the minimum energy (T100 (Fig- EI- The ratio cr/<7ioo 
is slightly lower than (6/n) 1 / 3 , this value corresponding 
to cubic precipitates showing only [100] facets. 

Robson and Prangnell§i deduced from experimental 
observations of A^Zr coarsening a Al/A^Zr interface 
free energy a = 100 mJ.m -2 at 773 K. The agree- 
ment between this value and the one deduced from our 
atomic model is perfect. In the same way, Hyland 22 
obtained from measured nucleation rates and incuba- 



tion times an experimental AI/AI3SC interface energy 
a = 94 ± 23 mJ.m" 2 between 563 and 623 K. Neverthe- 
less, this experimental value should be considered only as 
an order of magnitude as experimental nucleation rates 
and incubation times are hard to obtain. One has to be 
sure that precipitates of the critical size can be observed 
and the difference between the interface energies deduced 
from the incubation times or from the nucleation rates 
could be due to a detection limit for small precipitates 
greater than the critical size. Moreover, the Sc diffu- 
sion coefficient used by Hyland in his study differs from 
the one which has been more recently obtained from ra- 
dioactive tracer diffusion measurements^ and this would 
influence too the value of the interface energy deduced 
from his experimental observations. With these consider- 
ations in mind, this experimental value although slightly 
lower than the one we calculate (~113 mJ.m -2 ) is in 
good agreement with it. This indicates that the use of 
Wulff construction with mean-field theory is a good way 
to estimate this isotropic interface free energy and that 
our set of atomic parameters (Tab.QJ is realistic to model 
solid solution / Al 3 Zr and Al 3 Sc interfaces. 

C. Cluster size distribution 

For a dilute solution, the probability to observe in the 
solid solution a cluster containing n atoms having LI2 
structure ig^l 

C n ~ ^ = exp (-AGJkT), (11) 

where the formation energy AG n is given by Eq. JSJl. 
If the solution is supersaturated, the energy AG n is 
decreasing for sizes greater than the critical size and 
Eq. i|ll[l is assumed to be checked only for n < n* . As 
this is the criterion we chose to discriminate the solid so- 
lution from the LI2 precipitates (cf. section Hi C|l . this 
means that only the cluster size distribution in the solid 
solution should obey Eq. (|11J) and not the size distribu- 
tion of stable precipitates. 

We compare the cluster size distribution given by the 
Eq. I|llfl with the ones measured in Monte Carlo simula- 
tions for different temperatures between 723 and 873 K 
and different concentrations of the solid solution in the 
Al-Zr (Fig. |SJ as well as Al-Sc systems, both systems 
leading to the same conclusions. 

For stable solid solutions (s^ < x x)> au ener g e tic con- 
tributions entering AG n are positive and the cluster crit- 
ical size is not defined. Therefore, one expects Eq. (|llfl to 
be obeyed for all values of n. A comparison with Monte 
Carlo simulations shows a good agreement. The com- 
parison can only be made for small clusters: the proba- 
bility to observe large clusters in the simulations is too 
low to obtain statistical information on their distribution 
in a reasonable amount of computational time. This is 
interesting to note that for a solid solution having a con- 
centration equal to the solubility limit (x\ = x^), the 



8 



r 7r = 0.0274% 
.0548 % 

,2% 
5% 




u 



10"' 



10"" 



10" 



10"- 



10 H 





* = 2s 




t = 36 s 




f = 533s x 




f= 3290 s = 




f=2s 




f=36s 




* A A A a 








\ ( = 533 s x K 


r = 3290s X 






\m m X > i * f i 1*1 n n Si 



10 



»Zr 



"a 



FIG. 5: Dependence with the nominal concentrations x% r of 
the cluster size distributions of an aluminum solid solution 
at T = 773 K. At this temperature, the solubility limit is 
x^ t — 5.48 x 10 -4 . Lines correspond to prediction of clas- 
sical nucleation theory combined with CVM calculation and 
symbols to Monte Carlo simulations. 

prediction ^] of the cluster size distribution is still cor- 
rect. As the nucleation free energy is null for this con- 
centration, the only contribution to AG n arises from the 
interface. This shows that our estimation of the interface 
free energy a is coherent with its use in Eq. and that 
the capillary approximation gives a good description of 
the solid solution thermodynamics. 

For low supersaturated solid solution (for instance, on 
Fig. x\ x = 0.2 %), we observe a stationary state dur- 
ing Kinetic Monte Carlo simulations: the computational 
time to obtain a stable LI2 cluster is too high and the 
solid solution remains in its metastable state. Therefore, 
we can still measure the cluster size distribution during 
the simulations. The agreement with Eq. is still cor- 
rect (Fig.[3J). One should notice that now, the critical size 
being defined, the comparison is allowed only for n < n* . 

For higher supersaturations, the solid solution concen- 
tration xx is decreasing meanwhile stable precipitates 
appear (cf. kinetics of precipitation of AI3SC in Fig.[5J). 
This involves that the nucleation free energy is decreasing 
in absolute value and that the critical size n* is increas- 
ing. At each step we have to re-calculate the solid so- 
lution concentration and the critical size self-consistently 
by means of the definition Q of xx and by imposing that 
AG n (xx) is maximum in n* . Then we use this new value 
of the solid solution concentration in Eq. to calculate 
the corresponding cluster size distribution and compare it 
with the kinetic Monte Carlo simulation (cf. cluster size 
distributions in a Zr supersaturated aluminum solution 
in Fig. 0). We see that the time evolution of the cluster- 
size distribution is well reproduced by Eq. (|Tl"1) when the 
instantaneous concentration is used to calculate the nu- 
cleation free energy and therefore the prediction of cluster 
size distribution is not only verified during the nucleation 



FIG. 6: Evolution with time of the cluster size distribu- 
tions of an aluminum solid solution of nominal concentration 
x% r = 1 x 1CT 2 at T = 723 K. At this temperature, the sol- 
ubility limit is — 2.90 x 10~ 4 . Symbols correspond to 
Monte Carlo simulations and lines to prediction of classical 
nucleation theory combined with CVM calculation with the 
following instantaneous solid solution concentrations and crit- 
ical sizes: i Z r = lx 1CT 2 , 7x 1CT 3 , 2.7 x 10~ 3 , and 1.5 x 1CT 3 
and rz| r = 7, 8, 18, and 41 at respectively t — 2, 36, 533, and 
3290 s. 

stage but is well adapted even during the growing stage. 
Thus the thermodynamic description used in the classi- 
cal nucleation theory is in good agreement with results 
of atomic simulations. 



D. Kinetic description 

1. Diffusion 

Classical nucleation theory assumes that only 
monomers migrate and that larger clusters like dimers 
do not diffuse. We check this is the case with our atomic 
model by measuring during Monte Carlo simulations the 
diffusion coefficient associated to the gravity center of 
N atoms X in pure Al for 2 < N < 4. We obtain for 
Zr as well as Sc atoms that this diffusion coefficient is 
equal to the monomer diffusion coefficient divided by 
the number N of X atoms considered. This implies the 
following relation 

((f>x„) ^=f>r x „ 2 ), (12) 

where the brackets indicate a thermodynamic ensemble 
average and Arx„ is the displacement of the atom X„ 
during a given time. This relation is satisfied only if 
there is no correlation between the displacement of the 
N atoms X, which in other words means that cluster 
formed of the N atoms does not diffuse. In both sys- 
tems, tracer diffusion coefficient of Al is several order of 
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FIG. 7: Variation with nominal concentration and tempera- 
ture of the steady-state nucleation rate J st for AbjZr (top) 
and AI3SC (bottom) precipitation. Symbols correspond to 
Monte Carlo simulations and lines to classical nucleation the- 
ory combined with CVM calculation. 



magnitude larger than tracer diffusion coefficients of X. 
The relationship of Manning 54 shows that in that case 
the correlation factor /xx is almost equal to the tracer 
correlation factor fx which is equivalents^ to Eq. (|12|l . 
Thus the assumption used by classical nucleation theory 
of a diffusion controlled by monomers is checked for both 
Al-Zr and Al-Sc systems although interactions with va- 
cancies are different for these two binary systems. This is 
not the case for all systems: when vacancies are trapped 
inside precipitates or at the interface with the matrix, 
small clusters can migrateiSSiiLSi which affects kinetics 
of precipitation. 



2. Steady-state nucleation rate 

The steady-state nucleation rate is then predicted to 
be given by the equation^, 



where N s is the number of nucleation sites, i.e. the num- 
ber of lattice sites, AG* is the nucleation barrier and 
corresponds to the free energy of a precipitate of critical 
size n* as given by Eq. JSJl, Z is the Zeldovitch factor and 
describes size fluctuations of precipitates around n* , 



Z = 



(AG nuc Y 



2jr{a 2 a) 3 / 2 VkT 



(14) 



and (3* is the condensation rate for clusters of critical size 
n* . Assuming the limiting step of the adsorption is the 
long range diffusion of Zr or Sc in the solid solution and 
that Al atoms diffuse infinitely faster than Zr or Sc atom, 
the condensation rate is£ 



P* = -32tt- 



D 



x 
x x . 



(15) 



Although only monomers diffuse, the concentration ap- 
pearing in Eq. I|15|) is the nominal one as it reflects the 
gradient of concentration driving diffusion. Each time 
one Zr or Sc atom condensates on a cluster, three Al 
atoms condensate too on the same cluster. Thus clusters 
are growing from sizes An to 4(n +1). 

Comparing with the steady-state nucleation rate mea- 
sured in Monte Carlo simulations for different tempera- 
tures and different supersaturations of the solid solution 
in the Al-Zr and Al-Sc system, we see that the classical 
nucleation theory manages to predict J st (Fig. [TJ). The 
agreement is really good for low nominal concentrations 
of the solid solution (xx < 1 x 10~ 2 ) and is still good for 
higher concentrations: there is a small discrepancy but 
the relative values for different temperatures at a given 
concentration are correctly predicted. For instance, for 



the nominal concentration : 



Zr 



0.01, we obtain that the 



J st = N S Z{3* exp (— AG*/fcT), 



(13) 



steady-state nucleation rate is higher at T = 823 K than 
at T = 773 or 873 K. This shows that the kinetic model 
used by the classical nucleation theory is checked both 
for A^Zr kinetics of precipitation where there is repul- 
sion between the vacancy and the precipitating element 
and for AI3SC kinetics where there is attraction. 



IV. CAPILLARY APPROXIMATION 

Although it manages to catch thermodynamics of the 
solid solution, the capillary approximation that we used 
previously can look rough. First of all, one can wonder if 
it is reasonable to assume spherical precipitates especially 
for small ones. Moreover, when counting precipitates in 
Monte Carlo simulations, we assume them as being sto- 
ichiometric whereas in the mean-field calculation of the 
nucleation free energy we include anti-site defect contri- 
bution. Another source of mistake could be the use of the 
Wulff construction to calculate an isotropic interface free 
energy: doing so, we calculate the interface free energy 
of the most stable precipitate and therefore neglect some 
configurational entropy. 

In the present section, we calculate the cluster free 
energy without using the capillary approximation. The 
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results obtained with this direct calculation are then con- 
fronted to the ones obtained with the capillary approxi- 
mation. We also take the benefit of the exact results to 
discuss different levels of mean-field approximation used 
for the calculation of parameters entering in the capillary 
approximation. 



A. Direct calculation of cluster energies. 

Instead of using the capillary approximation to calcu- 
late the formation energy of Ll 2 clusters, we can calculate 
this quantity exactly. This can be done, following Ref.lo9l 
by sampling thermodynamic averages with Monte Carlo 
simulations so as to compute the free energy difference 
between a cluster of size n and one of size n + 1 at a 
given temperature. This method presents the drawback 
that a calculation is needed at every temperature of in- 
terest. We prefer calculating all coefficients entering the 
partition function as done in Ref. [13 and then derive the 
free energy at every temperature. 

A LI2 cluster containing nx X atoms can have dif- 
ferent shapes which we group by classes a of same en- 
ergy: D nXj0l is the number per lattice site of clusters 
containing rt x X atoms and having the energy H nXt0l — 



6e 



(i) 



6e 



(i) 



3c 



(2) 



3e 



.2) 



Energies are defined referred to the pure Al reference 
state as the cluster energy is the energy change due to 
the presence of a cluster in pure Al. The free energy of 
a LI2 cluster containing nx X atoms is then defined by 



G„ 



n x 



exp (—H nXta /kT) 

\ UU + b£ XX _ be AlAl + 6e XX ~ 6e Al 

fcTln ^£> nx , Q exp(-5F„ x , a /fcT)j (16) 



Degeneracies D nX;Q can be computed for a given size 
by generating clusters with a random configuration and 
then by counting for each energy level a the number of 
different clusters. The obtained values are presented in 
table Hill for LI2 clusters containing less than 9 X atoms. 
For bigger clusters, the degeneracy of the different classes 
is becoming too high to be countable. This is important 
to notice that we use the same criterion to define LI2 
clusters as in the kinetic Monte Carlo simulations (c/. 
section III Cfl , and that Zr or Sc atoms belonging to a 
LI2 cluster only have Al atoms as first nearest neighbors. 
We thus do not allow anti-site defects on the majority 
sub-lattice. This is not an important restriction as these 
defects have a high formation energy and therefore their 
contribution to the partition function can be neglected. 
As for the minority sub-lattice, anti-site defects can not 
be taken into account as they lead to a change of the 
precipitate size. 



TABLE III: Degeneracies D nx<a corresponding to classes 
of LI2 clusters containing nx X atoms and having en- 



ergy H nyL , a = n x [12uj w + 6e 
SH nXiCt for 1 < nx < 9. 
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The formation energies entering Eq. I|ll|l to calculate 
the cluster concentrations are the formation energies rel- 
ative to the solid solution, 



AG„ X (4) = G„ x - 2n xi u(x x ), 



(17) 



where fJ-(x x ) = (u x (x x ) — ^ai(^x)) is the effective 
potential, i.e. a Lagrange multiplier imposing that the 
nominal concentration of the solid solution is equal to 
the concentration of solute contained in the clusters as 
given by Eq. . As in the capillary approximation, this 
formation energy can be divided into a volume and an 
interface contribution: 



AG„ X (4) =4n x AG"" c (4) 

+ (367m x 2 ) 1 ^ 3 a 2 er Ilx , 

where we have defined the nucleation free energy 
~'(x x ) = (l2w« +64 



(18) 



AG" 



e (1) 
fc xx 



6c 



(i) 



+3e (2) 



3c 



(2) 
A1A1 



A1A1 

2u(x x )) /4 



(19) 



and the interface free energy 

-kT (36nn x 2 y 1/3 



a 2 cr ri 



In 



Dn x ,a exp (~~8H nXta /kT) 



(20) 



All information concerning the solid solution, i.e. its 
nominal concentration, is contained in the nucleation free 
energy whereas the interface free energy is an intrinsic 
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FIG. 8: Variation with the nominal concentration x\ Y of the 
nucleation free energy AG nuc at T = 723 K obtained with dif- 
ferent approximation: direct calculation of the cluster forma- 
tion free energy (Eq. 1191 ) or capillary approximation with the 
nucleation free energy given by the CVM calculation, the ideal 
solid solution model, and the regular solid solution model. 



property of clusters, which was already the case with the 
capillary approximation. The main difference is that now 
the interface free energy depends on cluster size. Perini 
et alw%> show that this size dependence can be taken into 
account in the capillary approximation by adding terms 
to the series (JSJ of the formation energy reflecting line 
and point contributions. 

We compare the nucleation free energy obtained from 
this direct calculation of the cluster formation energies 
with the one that we previously calculated with CVM 
in section IIII Al (Fig. ISJ . The direct calculation leads to 
a slightly lower nucleation free energy in absolute value 
than the CVM one. This mainly arises from the neglect 
of excluded volume between the different clusters in the 
direct calculation. Nevertheless, the agreement is cor- 
rect for all temperatures and for both Al-Zr and Al-Sc 
systems. This shows that these two approaches used to 
describe thermodynamics of the solid solution, i.e. the 
mean-field and the cluster descriptions, are consistent. 

The interface free energy defined by Eq. I|2U|I is de- 
creasing with cluster size, the variation becoming more 
important at higher temperatures (Fig.|5J). The asymp- 
totic limit is smaller than the interface free energy a that 
we calculated in section IIII Bl using Wulff construction 
and Bragg- Williams approximation. This is quite natural 
as Wulff construction predicts the cluster shape costing 
least energy. We are thus missing some configurational 
entropy by using it to compute an interface free energy 
a and we overestimate a. This error can be neglected at 
low temperature (T < 773 K) where precipitates show 
sharp interfaces but it increases with temperature when 
precipitate shapes are becoming smoother. 

We use this direct calculation of the cluster formation 
free energy (Eq. (|16|) and l |17f l) to predict cluster size 
distributions in the solid solution and compare the re- 



FIG. 9: Variation with the cluster size nzr of the interface 
free energy between the solid solution and AFjZr. Symbols 
correspond to a nZr as given by the direct calculations of the 
cluster formation free energy (Eq. 1201 ) and lines to a, i.e. 
the Bragg- Williams calculation combined with the Wulff con- 
struction. 
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FIG. 10: Cluster size distribution of two aluminum solid solu- 
tions of nominal concentrations x| r = 8 x 10~ 4 and 2.4 x 10~ 3 
at T = 873 K. At this temperature, the solubility limit is 
= 1.6 x 10~ 3 . Symbols correspond to Monte Carlo simu- 
lations and lines to prediction of classical nucleation theory as 
given by Eq. 11 U . To evaluate the cluster free energy of for- 
mation, we use the capillary approximation (Eq. (JHJ) with the 
nucleation free energy given by CVM for the continuous line 
and the direct calculation (Eq. I|16|l and $TQ ) for the dashed 
line. 

suits with the distributions obtained with the capillary 
approximation (Eq. JHJ) combined to the CVM calcu- 
lation. These two models lead to similar distributions 
fFig. 1 1 Of) . indicating that the associated thermodynamic 
descriptions are consistent. Nevertheless, the distribu- 
tion predicted by the direct calculation better reproduces 
the ones measured during the Monte Carlo simulations. 
Thus, the capillary model is good to describe thermody- 
namics of the solid solution but it can be improved. 

Comparing the steady-state nucleation rates predicted 
by the two thermodynamic models with the ones mea- 
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FIG. 11: Variation with nominal concentration x^ r and tem- 
perature of the steady-state nucleate rate J st for AlsZr. Solid 
lines correspond to prediction of the classical nucleation the- 
ory when using the capillary approximation with the CVM 
calculation of the nucleation free energy and dashed lines 
when using the direct calculation of the cluster formation free 
energy. Symbols are measurements in Monte Carlo simula- 
tions. The errorbars correspond to the uncertainty on the 
measurements of J st due to the choice of the critical size cor- 
responding to each energetic model. 



sured in Monte Carlo simulations (Fig. 111( 1. we do not 
obtain any improvement by using the direct calculation 
of cluster energy instead of the capillary approximation. 
For low supersaturations, both models are in reason- 
able agreement with Monte Carlo simulations whereas for 
higher supersaturations discrepancies appear. The direct 
calculation leads to a slightly lower nucleation rate than 
the capillary approximation. This mainly arises from a 
difference of the critical size: is usually 1 atom greater 
with the direct calculation than with the capillary ap- 
proximation. As the use of the direct calculation im- 
proves the agreement for the cluster size distribution, the 
discrepancy observed at high supersaturations is not due 
to a bad description of the solid solution thermodynamics 
but may arise from limitations of classical nucleation the- 
ory itself. The assumption of a constant flux between the 
different size classes made by this theory to solve the rate 
equations associated to the cluster size evolution may not 
apply at high supersaturations. This can be seen in our 
atomic simulations by the fact that, for these supersatu- 
rations, the linear domain observed for the variation with 
time of the number of precipitates and used to define the 
steady-state nucleation rate is more restricted than in 
the low supersaturation case shown on Fig. ^ One could 
try to improve the agreement with atomic simulations by 
using more sophisticated mesoscopic models like cluster 
dynamics£2*£ii§2iS which do not need such a kinetic as- 
sumption to solve the rate equations. Another improve- 
ment that could be made to classical nucleation theory is 
to consider the variation with the nominal concentration 
of the diffusion coefficient of X atoms which would lead 
to a diffusion coefficient different from the impurity one 



that we use. 



B. Other mean-field approximation 

Usually, one does not calculate the nucleation free en- 
ergy with CVM as we did in section ITlI Al but one uses 
simpler mean-field approximation to evaluate the chem- 
ical potentials entering Eq. © of AG nuc . We test these 
other approximations and see if they are reliable to be 
used with classical nucleation theory. 

The easiest approximation that can be used is the ideal 
solid solution model in which one keeps only the configu- 
rational entropy contribution in the expression of chemi- 
cal potentials and calculates this term within the Bragg- 
Williams approximation. This leads to the following ex- 
pression 



MB! 



i 



In 



(21) 



The exact expression of the nucleation free energy, i.e. 
with the cnthalpic contribution, can be calculated within 
the Bragg- Williams approximation too. This is called the 
regular solid solution model and gives 



AG£^(4) 



a smuc 
^^idcal 



(X 



eq2 



+ i ((i _ 4«) 2 - (i - 4) 2 ) 



(22) 



Comparing all different mean-field approximations 
used to evaluate the nucleation free energy (Fig. |8J, we 
see that for low supersaturations all approximations are 
close, but that for an increasing nominal concentration of 
the solid solution discrepancies between the different ap- 
proximations are becoming more important. Both ideal 
and regular solid solution models overestimate the nu- 
cleation free energy compared to the CVM and the di- 
rect calculations. The discrepancy is even worse when all 
contributions, i.e. the enthalpic and entropic ones, are 
considered in the Bragg- Williams approximation. Thus, 
when the supersaturations is becoming too important, 
the Bragg- Williams approximation seems too rough to 
give a reliable approximation of the nucleation free en- 
ergy. 

This becomes clear when combining these approxima- 
tions of AG nuc with classical nucleation theory to pre- 
dict cluster size distributions. The ideal and the regu- 
lar solid solution models completely fail for high super- 
saturations to predict the cluster size distributions ob- 
served during Monte Carlo simulations (Fig. I12fl . The 
predicted critical size is too small as it corresponds 
to a cluster size in the observed stationary distribution 
and the predicted probabilities for each cluster size are 
too high compared to the observed ones. As the pre- 
diction of the steady-state nucleation rate by classical 
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FIG. 12: Cluster size distribution of an aluminum solid solu- 
tion of nominal concentration x$ c = 7.5 x 1CP 3 at T = 773 K. 
Symbols correspond to Monte Carlo simulations and lines 
to prediction of classical nucleation theory with the different 
mean field approximations of the nucleation free energy. 



nucleation theory is based on the predicted size distri- 
bution, the ideal solution model and the Bragg- Williams 
approximation lead to an overestimation of J st too. Thus 
the use of CVM to calculate nucleation free energy re- 
ally improves agreement with atomic simulations com- 
pared to more conventional mean-field approximations. 
This arises from the fact that order effects are not taken 
into account in Bragg- Williams approximation whereas 
they are in CVM. These order effects correspond to a 
strong attraction for first nearest neighbors and a strong 
repulsion for second nearest neighbors between Al and 
Zr atoms as well as Al and Sc atoms. They are the rea- 
son at the atomic scale why a supersaturated aluminum 
solid solution evolves to lead to the precipitation of a Ll 2 
compound. Therefore one must fully consider these order 
effects when modelling kinetics of precipitation. 



V. CONCLUSIONS 

We built an atomic kinetic model for Al-Zr and Al-Sc 
binary systems so as to be as close as possible to the real 
systems. Thanks to this model, we were able to simulate 
at an atomic scale kinetics of precipitation of the LI2 
ordered compounds A^Zr and AI3SC. 

From this atomic model we deduced the corresponding 
interface and nucleation free energies which, with the dif- 
fusion coefficients, are the only parameters required by 
mesoscopic models like classical nucleation theory. When 
CVM is used to calculate the nucleation free energy we 
showed that the capillary approximation leads to a satis- 
fying thermodynamic description of the solid solution. If 
one wants to improve this description, one can calculate 
directly formation free energies of the different size clus- 
ters. This leads to a better description of the thermody- 
namic behavior of the solid solution, as the agreement on 



cluster size distribution is better, but it does not dramati- 
cally change predictions of the classical nucleation theory. 
This shows that the capillary approximation is reason- 
able. From the kinetic point of view, classical nucleation 
theory assumes that evolution of the different clusters is 
governed by the long range diffusion of monomers. For 
Al-Zr and Al-Sc systems, it appears to be a good as- 
sumption as we checked that di-, tri-, and 4-mers do not 
diffuse and that the steady-state nucleation rates mea- 
sured in Monte Carlo simulations are in good agreement 
with predictions of the classical nucleation theory. Dis- 
crepancies appear at higher supersaturations which may 
be due to the dependence of the diffusion coefficient with 
the solute concentration of the metastable solid solution 
or to the limits of the classical nucleation theory which 
requires the nucleation regime to be separated from the 
growth regime. Nevertheless, the nucleation model was 
built on purpose to predict kinetics at low supersatura- 
tions for which kinetic Monte Carlo simulations are not 
tractable. 

On the other hand, when one uses less sophisticated 
mean field approximation than CVM like the Bragg- 
Williams approximation to calculate the nucleation free 
energy, predictions of the classical nucleation theory com- 
pletely disagrees with Monte Carlo simulations, espe- 
cially when supersaturations are too high. This shows 
that short range order effects which are naturally consid- 
ered in CVM must be taken into account so as to build 
a kinetic mesoscopic model based on a reasonable phys- 
ical description. This is expected to be the case for all 
systems where order effects are important and thus for 
systems leading to the precipitation of an ordered com- 
pound. 
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APPENDIX: WULFF CONSTRUCTION 

We use the Wulff constructionist so as to define an 
isotropic interface free energies a from the free energies 
cioo, cr iio ; and (Tin- This construction allows us to de- 
termine precipitate real shape and to associate with it 
a which corresponds to the same interface energy for a 
spherical precipitate having the same volume. 

AI3X precipitates will show facets in the [100], [110], 
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and [111] directions if the following conditions are met: Considering a spherical precipitate with the same volume 

and the same interface energy, one gets 

\/2/2 ctioo < ciio < cioo, 

2^6/3 cuo — V3/3 ClOO- 

For Al 3 Zr and Al 3 Sc, with the set of parameters given 
by table HI this is true for all temperatures. Each facet 
surface will then be proportional to 



Tioo = 4((7ioo - \/2cno) 2 c = y ^-(6cioorioo + 12cii rno + 8c m r m ). (A. 4) 

-2(cioo -2x/2a 110 + %/3c m ) 2 , (A.l) 
Tno = 2V2(-2cioo + V2(Tii ) 

{V2<r 110 - V3c m ), (A.2) 
r m = 3V3/2 (-Ci 00 _ 2o no + a ui) 

+3/2CT 10 o(4\/6cii -6c m ). (A.3) 
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